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The  purpose  of  this  paper  is  to  present  an  overview  of  several 
n-body  dynamics  formulations  in  the  spacecraft  dynamics  literature. 

Even  though  the  emphasis  is  different,  the  overview  in  this  paper  is 
somewhat  in  the  spirit  of  Likins  (1970,  1974,  1975)  and  Meirovitch  (1975), 
and  the  background  references  for  this  paper  are  essentially  the  same 
as  for  these  earlier  overviews.  This  paper  differs  from  previous  papers 
in  the  spacecraft  dynamics  literature  in  that  both  ''momentum  formula- 
tions" and  "velocity  formulations"  are  discussed  in  a single  language- - 
the  language  of  the  transformation  operator  formalism  (Jerkovsky,  1976). 

The  paper  starts  out  with  a description  of  multi-body  tree  config- 
urations. A path  matrix,  n,  and  a reference  matrix,  c , are  defined  in 
the  spirit  of  Roberson  and  Wittenburg  (1966).  The  next  step  is  the  intro- 
duction of  the  "primitive"  or  "free  body"  equations  of  motion  in  terms 
of  a single  equation.  This  equation  is  then  linearly  transformed  via  a 
transformation  operator.  A,  and  the  result  is  a new  "transformed" 
equation  of  motion.  This  method  of  transforming  "old"  differential 
equations  to  "new"  differential  equations  is  based  on  Kron's  method  of 
subspaces  (Hoffmann, 1944),  and  is  similar  to  the  matrix  method  of 
structural  analysis  (Przemieniecki,  1968).  In  the  old  differential  equa- 
tions the  velocities  are  inertial  velocities,  whereas  in  the  new  differen- 
tial equations  the  velocities  are  relative  velocities.  The  transformation 
to  relative  velocities  is  made  so  that  relative  velocity  constraints  can  be 
treated  more  readily;  this  transformation  is  made  in  the  spirit  of  class- 
ical mechanics  (Corben  and  Stehle,  1960)  where  generalized  coordinates 
are  introduced  so  that  the  constraints  become  trivial. 

As  .an  alternative  to  transforming  to  relative  velocities,  the 
equations  of  motion  can  be  kept  in  terms  of  inertial  velocities,  and 
the  relative  velocity  constraints  can  then  be  incorporated  via  Langrange 
multipliers.  This  alternative  approach  is  particularly  attractive  in 
cases  where  it  is  not  a simple  matter  to  express  the  inertial  velocities 
of  a multi-body  system  in  terms  of  an  independent  set  of  relative  veloc- 
ities; such  a case  occurs  when  the  multi-body  configuration  is  not  a 
tree--i.e.  , when  there  .are  closed  loops. 
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'I'iic  I ii-His  j)  ri- s (■  I lU’ti  in  thi.s  p.ipii'  .iMsun.i-  that  tiic  mulli-body 

I onl  1 all  ra  t i on  lonsi  st  s of  n ripid  bodios.  '1  lio  -sann'  proc  i-dun-  i an  bo 

nsod  if  soir.o  or  all  tlu'  b>)dii  s a ro  lloxibli. i.o.,  tiio  .-,1  i' net  u ro  ol  iho 

muili-bod\  dynaniios  o{;nalions  is  Iho  same  \\hothor  tho  bodio.-^  ari- 
flosiblo  or  i-ioid.  In  fait,  a largo  nundjor  of  ri_;id  bodios  i an  bo  usod 
to  modol  a floxiblo  bofly,  and  Iho  striutnro  of  tiio  oi|ualions  duos  not 
dopond  on  I'.hothor  n is  hii'm-.  I;  it  is  rlosirod  to  troat  all  n borlios  as 
iloxiblo,  tluMi  thoro  aro  throo  modific.itions  v\hi<h  ,iro  ro(|iiirod:  (1)  tho 

'nrinutiv.o"  or  ''fri'o  body'  iMjU.ition  of  motion  must  im.!u<io  i.-quations 
ol  motion  for  the  doform.ition  dooroos  ol  Iroodom  (Hodloy  ot  al.  , l‘>72; 

Jorkov  sky,  l9T7al;  (2)  tin-  "primitivo"  and  ''t  ransformi  d'  \ olooilios 
nui.-^t  iiu  ludo  tho  limo  dorivati\  os  of  tho  doformation  i oordinati's 
(liodtoy  ot  <al.  , Ih75);  .md  (’<)  tho  transformation  oporator,  A,  must 
oxpross  tho  inortial  \olooitios  in  torms  ol  rolativo  \oloi  itios  ])lu.-' 
di'lorm.ition  t oordinati’s  timo  dorivativos. 

Tho  nioniontum  formulation  and  vadooity  foririulation  oqualions 
dosoribod  horoin  aro  similar  to  the;  equations  that  might  bi>  obtainod 
using  a I l.imillonian  or  L.igrangian  mechanics  approach,  rospocti  voly . 
Ilowovor,  there  are  t\\.o  fundamontal  differences:  (1)  In  tho  present 

a))proach  tho  Hamiltonian  (a  function  of  generalized  loordinatos  plus 
gone  r.i  lizi'd  momenta)  or  Lagrangian  (a  function  of  generalized  coordi- 
nates plus  tho  time  derivatives  of  these  generalized  coordinates)  aro  not 
■ormi'd,  and  the  equations  of  motion  are  not  obtained  by  ]>artial  differ- 
r.tiation  of  tho  Hamiltonian  or  Lagrangian:  (2)  In  the  pre-sent  approach 
I'in.il  ocpiations  are  not  oxprc'ssed  c'X]ilicitly  and  solely  in  terms  of 
gciior.ilizod  coordinates  and  gonor.ilizod  mcinumta  or  tho  time  deriva- 
tives of  tho  goiiorali/.od  coordin.atc's;  instinid,  the  final  oc|u;iticins  are 
expressed  in  terms  of  some  "i  nt  o r mediat  e"  or  "auxiliary"  variables 
which  ,1  ro  algebraic  functions  of  tho  generalized  coordinates  and 
gone  ralizofj  momenta  or  the  time  derivatives  of  the  generalized  coordi- 
nati's.  Iho  intermediate'  or  auxiliary  variables  generally  have  physical 
sigiulicatK  o.  In  tlu'  i asi'  of  the  momc-ntum  formulation,  the  intermc'di- 
ato  variables  aro  the  velocities  of  hinge  points  and  the  translation.il 
momenta,  even  if  there  ,ire  no  relative-  translational  degrees  of  freedom. 
The  justifii  ation  for  the  retention  of  intermediate  variables  is  the  con- 
ceptual .end  c omputational  sinqilicity  ot  the  resulting  equations.  Note 
that  a Hamiltonian  formulation  vvcjuld  not  even  include  any  time 
derivative's  of  gene  ralized  coordinates  (the  generalized  momenta  would 
be-  used  instc-ad).  Hence,  a momentum  formulation  which  retains 
velocities  is  re-ally  a "mixed"  formulation  rather  than  a "pure"  formu- 
lation. But,  similarly,  a vc'locity  formulation  which  retains  velocities 
otlter  than  the  time  derivatives  of  the-  gcneralizc-d  coordinates  is  also 
a n ixed  formul  ition.  .Note  that  if  a pure  formulation  is  used  then  the 
terms  in  the  i.cpiations  of  motion  are  unique  (if  properly  symmetrized); 
on  the  othc'i-  hand,  in  a mixed  formulation  the  terms  depend  on  the 
particular  choice  "f  intermediate-  vnriablc-s. 
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Dt  sf  ription  Of  Mulli-Bod y Tree  Coiifigurationei 

Givon  an  n-body  configuration,  we  label  tlie  bodies  from  1 to  n, 
assigning  the  label  1 to  the  "rnain''  or  "central"  body.  Wc  obtain  a 
"graph"  of  the  configuration  by  putting  each  body  in  correspondence  with 
a verte^x  (or  node)  of  a graph  and  connecting  any  two  vertices  of  this  graph 
with  a branch  if  the  corresponding  bodies  have  any  degrees  of  relative 
motion  between  them.  If  the  resulting  graph  is  a tree  (i.e.  , if  there  are 
no  closed  loops),  then  the  n-body  system  is  said  to  have  a tree  configura- 
tion . If  the  graph  is  not  a tree,  then  a tree  can  still  be  associated  with 
the  graph  by  cutting  as  many  branches  as  there  are  closed  loops.  Any 
branch  in  any  closed  loop  may  be  cut,  and  different  choices  will  lead  to 
different  trees  of  the  graph. 

Thus,  to  any  n-body  configuration  there  corresponds  a tree  with 
Body  1 at  the  center  of  the  tree.  For  the  moment  we  do  not  concern 
ourselves  with  how  many  degrees  of  freedom  there  are  between  adjacent 
vertices  (i.e.,  between  adjacent  bodies),  or  if  the  actual  n-body  configu- 
ration has  closed  loops  or  not. 

We  label  the  bodies  (or  vertices  of  the  tree)  such  that  all  the  bodies 
between  Body  1 and  Body  j have  an  index  i between  1 and  j.  Also,  we 
let  j be  the  set  of  integers  which  includes  1 and  j and  also  includes  the 
labels  of  all  bodies  between  Body  1 and  Body  j.  Let  be  the  label  of  the 
body  next  to  Body  j on  the  path  from  Body  j to  Body  1;  similarly,  let 
be  thelabelof  the  body  next  to  Body  j;  etc.  Then,  the  set  j"  consistof 

the  labels  j^=  (j,  j,  j l}.  Evidently,  ^ is  the  set  of  indices  of  all 

bodies  which  are  "inward"  from  Body  j,  including  Body  j.  An  example 
9-l^ody  tree  configuration  is  shown  in  Fig.  1.  Fig.  2 shows  the  sets 
to  9 for  this  example.  The  labels  i to  9 are  also  shown,  where  has  beer 
set  to  zero  so  that  j is  defined  for  j = l to  n. 

Next,  we  let  k^be  the  set  of  all  j such  that  k is  contained  in  j 
(i.  e,  , such  that  k C j).  Evidently,  k is  the  set  of  indices  of  all  bodies 
which  are  "outward"  from  Body  k,  including  Body  k,  _Fig.  2 shows  the 
sets  1 to  9 for  the  9-body  example.  Note  that  the  set  1 includes  the  integers 
from  1 to  n ( 1 = { 1,  2,  . . . , n })  because  all  bodies  are  outward  from 
Body  1. 

We  now  introduce  the  "path  matrix"  as  follows.  Letting 
denote  the  element  of  n in  the  i row  and  j^^  column,  we  define 

1,  if  Bodyj  is  between  Body  1 and  Body  i 
0,  otherwise 

Fig.  2 shows  for  the  9-body  example.  Note  that  the  rows  of  ^ are 
determined  by  the  inward"  sets  i for  i 1 to  9,  and  the  columns  of 
•■are  determined  by  the  "outward"  sets  j . 
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If  we  think  of  Body  j as  the  body  to  which  Body  j is"referenced" , 
then  we  can  doline  the  "reference  matrix"  C'as  follows 


r . . 
1.1 


1>  if  i ■=  j 

-1,  if  Body  i is  referenced  to  Body  j 
0,  otherwise 


(2) 


Fig.  2 shows  ‘ for  the  9-body  example.  Evidently,  the  -1  terms  of 
can  bo  filled  in  by  looking  at  the  labels  2.  3,  ....  n.  Note  that  the 
labels  2,  3,  ....  n define  the  tree  configuration  completely,  because 
these  labels  allow  us  to  draw  the  tree  and  then  the  sets  j and  k for 
j,  k = 1,  2,  ....  n.  Also  note  that  column  j of  p shows  which  bodies  are 

referenced  to  Body  j. 

From  Fig.  2 we  note  that  nand  pare  lower  triangular,  and  p 
is  sparser  than  n However,  rr  and  p contain  the  same  information,  in 
fact,  Ti  and  rare  inverses  of  each  other: 

pn=f=np  (^) 

n 

where  1 is  the  nXn  identity  matrix.  We  will  see  later  that  after  the 
matrice’s  ~ or  P are  introduced  for  a particular  tree  configuration, 
the  equations  become  independent  of  the  particular  configuration  under 
consideration. 

We  now  define  the  set  i j to  be  the  intersection  of  the  sets  i and  j ; 


i j = i n j = j i 


(4) 


Now  all  the  sets  k,  for  k = 1 to  n,  are  "nested"  in  the  sense  that  if  the 
sets  I and  J have  any  elements  in  common,  then  either  i is  contained  in 
J,  or  j is  contained  in  i.  Hence,  i j is  either  i or  j or  the  empty  set  0. 
Symbolically , 


i j 


i , if  i C j 
j . if  i C i 
0,  otherwise 


(5) 


Note  that  i i = i. 

We  now  introduce  the  points  ci,  hj,  and  bj  as  follows,  for  i = 1 to  n 
Let  Ci  be  the  center  of  mass  of  Body  i.  If  Body  i is  a rigid  body,  then 
the  point  ci  is  a fixed  material  point  of  Body  i;  if  Body  i is  deformable, 
then  c "floats"  in  the  body.  Whether  Body  i is  rigid  or  deformable, 
i 
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lot  li.  bi‘  a fixed  material  f)oint  which  is  the  "hinge"  point  for  Body  i. 

Let  ij.  be  a fixed  material  point  in  iiody  i to  which  Body  i is  referenced; 
the  point  b^  is  the  "Ijase"  point  for  Body  i.  By  convention  bj  is  tlie 
inertial  reference-  origin  (i.c-.  , a point  fixed  in  "liody  0").  Fig.  1 shows 
the  points  cj,  hj,  and  bj  for  the-  V-body  example. 

Primitive  Figuations  of  Motion 


By  "[irimitive"  eciuations  and  variables  we  will  mean  equations 
and  variables  which  refc-r  to  each  body  as  a separate  and  distinct  body, 
without  regard  to  how  it  fits  into  the  multi-body  configuration. 

Let  P^  be  the  liru-ar  momc-ntum  of  Body  i,  and  let  !■' ^ be  the 
force  on  Boc^  i.  Let  1 h*  j be  the  <ingular  momi-ntum  ol  BocJy  i about 
c^,  ana  let  Lj-  be  the  torque  on  Body  ' about  cj.  Note  that  we  arc'  using 
the-  body  index  i as  a supersc  ript,  and  the  point  i j as  a sulisc  ript.  Lc-t 
Vej  be-  tile  linear  velocity  of  the  point  cq;  thus,  if  q-  is  th<-_  position 
vec  tor  to  c^  from  the-  inertial  referc'uce  origin,  then  v^  . “ ^ whc're 

the  dot  ovc-r  a vector  is  used  to  denote  the-  time  derivative  in  the  inc-rtial 
retorc-nce  frame-.  L,c-t  ^ be  tlu-  anguhir  velocity  ol  a frame  fixed  in 
Body  i.  For  simplic  ity,  we  now  assume  th.at  all  bodic-s  are  rigid;  we 
have-  alrc-ady  i- lent  ioiu-d  t he  modifi  c a t ions  recjuin-cl  whe-n  .some-  or  all 
bodies  are  deformable.  Let  M*  be  the  mass  of  Bc.)dy  i,  and  let^^-  be 
Uu-  inerMia  (dyadic)  of  Body  i about  cj;  let  VV*^  = (M^)*l  and 

Now  dc-finc-  G,  K,  , and  Y to  be-  column  matrices  of  2n 
Gibbsian  ca-clors  as  follows  (Jc-rkovsky,  1976) 


wlu-n-  0 isj.hc-  ve^c  tor,  and  FiJ-  is  the-  "Euler  coupling  force" 

on  ihiciy  i:  E‘.  ■-  -P  • ll^j  . Also,  define  (J.  and  - as  diagonal  matrices 

of  positi\e  cic-finite  symmetric  dyadic  s as  follows 


n -1 

I 

Y2 

> 1 

kV 

w' 

w? 

n 

n 

M 

W 

whore  1,  v\'^  = VV^  I and  1 is  the  idenlit^^dyadie . The  off - 

diagonal  (.•leinents  v>l  -^and  • are  thi-  zero  dyadic  O . The  primitive 
momentum  formulation  ecpiations  for  the  n-body  system  are  now 
;;i\en  1)^ 

G X - K.  Ci  = tJ.  • or  = ■/  • G (8) 

wlu  re  X IS  .1  lolumn  matrix  of  /.ito  vectors  which  is  included  here 
only  for  |u-dacoai<-  il  reasons.  G,  K,  aiifl  will  be  called  the  primitive 
system  niomentum,  lone,  and  \ ('locity,  respi'i  tively . U and  will 
be  tailed  the  primitive  system  mass  and  inverse  mass,  respectively. 
The  \a  lot  ity  formulation  ecjuation  is 


-1  • r V = K (^) 

The  kinetii  enert’y  of  the  system  t>f  n rigid  bodies  is  given  by 

T = i c/  ■ ~ ^ • C (in) 

where  G^  and  ' are  the  transijoses  of  G and  ; thus, 

G^^  .lV  ll''^  f'  ■■■  P'l]  yirnilarly  for  0-^  The 

time  deri\ativ(‘  ol  the  kinetii"  energy  is  given  by 

T = k‘  • ■ (11) 

In  the  next  set  tion  we  express  the  primitive  system  velocity  " in  terms 
of  a new  (or  transformed)  system  velocity  ' . This  then  automatically 
defines  new  system  variables  so  that  Eqs.  (8)  to  (11)  maintain  their  form. 

Transformed  Equations  of  Motion 


We  will  ultimately  be  interested  in  the  case  where'  there  are  less 
than  6 free  degrees  of  freedom  between  some  or  all  of  the  adjacent 
bodies  of  the  tree  configuration.  Therefore,  we  will  now  transform  to 
relative  velocities  which  can  then  be  prescribed  if  the  corresponding 
degrees  of  freedom  are  constrained. 
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Velocity  Transformation 

We  will  now  express.tht;  inertial  velocities  ^ and  v^.  in  terms 
of  the  relative  velocities  and  U^.  We  define  as  the  angular 
velocity  of  Body  i with  respect  to  Body  i.  to  which  it  is  referenced: 

iV  = - 7-*“  (12) 

0 ^ 

’’BodyO''  is  the  inertial  reference  frame,  and  hence  a = 0.  Therefore, 

Lc-t  - r^  - rjjjBe  the  position  vector  to  point  a from  point  b. 

Also,  for  k = 1 to  n,  let  v denote  the  time  derivative  of  the  vector 
V with  respect  to  Body  k;  then,  y - . Now  define  as  the  time 

derivative  with  respect  to  Body  j.  of  the  position  vector  to  the  hinge  point 
h^  in  Body  i from  the  base  point  bi  in  Body  i : 

i ~ 


Note  that  U - R 


, . = V,  since  r , 

hjbi  bj 


The  primitive  inertial  velocities  and  v^-.  can  be  expressed 

in  terms  of  the  transformed  relative  velocities  as  I'ollows  (Jerkovsky, 
1977  b) 


^ n.  V 

j = l 1 J J e a 1 J 


+ u-*) 


where  we  make  use  of  the  notation  that  for  any  two  vectors  and  P . 

^he  dyadic  of  O'  is  denoted  by  5 and  is  defined  by  S • o^=  5xp  , and 

is  the  dyadic  transpose  of  Ot  , and  is  defined  by  3:  • b = P • o*  ; note 
that  j:  t - a,  i.  e.  , is  skew-symmetric.  The  inverse  of  Eqs.(i4)  are 


= 2: 


TJJ  = Ir  ^ . iTi 


U = E f • • (R  u • '»■’  a V ) = R , . 

. , 11  c.h.  c.  c.h. 

J=1  J 1 J 11 

Note  that  Eq.  (15)i  is  the  same  as  Eq.  (12). 


av  ) = R , .a^  + v -R  ,-ai-v 

c.  c.h.  c.  c.h.  c. 

J 11  111  1 


We  now  define  * to  be  a column  matrix  of  vectors  (and,  therefore, 
is  a row  matrix  of  vectors)  as  follows: 


~=  B ■ c , 


where  B = A 


-1 


^ A • 


(17) 


A is  the  iranslormatiori  operator  which  expresses  the  primitive  system 
velocity  in  terms  of  the  transformed  system  velocity  ” , and  B is 
its  inverse.  A and  B are  matrices  of  dyadics;  their  elements  can  be 
obtained  by  inspection  of  Eqs.  (14)  and  (15).  Note  that  A • B = B • A 
is  an  identity  with  Zn  1 's  on  the  diagonal.  See  Jerkovsky  (1976,  1977b) 
for  examples  of  A and  B matrices. 


Induced  Transformations 


In  order  for  Eqs.  (10)  and  (11)  to  maintain  their  form,  we  define 
the  transformed  momentum,  force,  and  mass  as  follows: 


G = A • G , 


K = A 


K , 


p = A ■ p • A 


(18) 


Carrying  out  the  operation  A • G,  we  find  G as  follows: 


C ' = L h) 


where 


h: 


pi  p2 


(19) 


II 

E 

i = l 

(H^ 

1 1 c. 

+ R u 
c.h. 

J 

1 

1 J 

— 

n 

pJ  = 

E 

i = l 

- . . P^  = 
ij 

E p' 

ifj 

p')  = E.  (»c  ^ ^c  h ■ 

ifj  i i j 


(20) 


Evidently,  P'1  is  the  linear  momentum  of  the  set  of  bodies  whose  labels 
are  in  the  set  j ; i.  e.  , the  set  of  bodies  outward  from  Body  j,  j^cluding 
Body  j.  If  we  refer  tq_this  set  of  bodiejs  as  "System  j ",  then  PJ  is  the 
momentum_of  System  j.  Similarly,  H is  the  angular  momentum  of 


of  System  j about  the  hinge  point  hj  in  Body  j.  Carrying  out  the  operation 
a1  • K yields  the  same  type  of  equations,  with  L'i 
and  P's,  respectively. 


's  and  F's  replacing  H's 


The  definitions  of  G and  p given  in  Eq.  (18)  now  yield  G = p ■ a 
in  the  form 


H 


H 


f 1 . . 

. ^In 

^h,h 
1 n 

Cijhi 

•rhn 

h 

n n 

^c-jh 
nl  n 

s^  . • • 

"il^i 

■ h 

c-  ih 
nl  n 

nil 

M 

s* 

Cr-h, 

In  1 

■ 

nn  n 

S 


Cr-h 

In  1 


c-  -h 
nn  n 

M 


M 


nn 


1 


Q 


u‘ 


u 


n 


(21) 
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w tu-  ri.‘ 


T ■'  - 

h.h 
1 .1 

11 

i:  ■ 

k=i 

n k, 

* k 

. R , 

1\  1 

• r’ 

‘ k'\i 

11 

L ^ 

M R , • 
c , h. 
k 1 

K 

c.  h 
k J 

n 

] 

L 

'^c.  ,h. 
hi  t 

- NpJ 

It 

ki  " k, 

K J - 

‘k''i 

« H,. 

k < 1 1 

l^ 

(22) 

n 

i 

M ^ 

- "i  • 1-  1 

m‘'T  x:._ 

M"  *P 

k f i J 

Ri'vall  that  tlu’  set  i j is  the  inlersei  turn  of  the^_set,s  i ;i  n<l  J.  11  i j = 0, 

tluMi  thf  sum  over  k t ij  yields  tlu‘  zero  dyadic  O.  ll_i  j,  then 

It;' h I '*^'h|  > which  is  the  inertia  (dyatiii  ) of  System  j about  hj;  ,ilso, 

M.ij  MJ , which  is  the  mass  of  Systimi  j.  Note  that  Rijjhj 

position  vei  tor  to  the  (loint  >d  j , the  center  ol  mass  oi  Systion  i j , irom 

till'  point  hj . For  an  I'sampK'  ol  , see  Jerkovsky  (ld7(),  l‘t7ib). 


Ecjs,(18l  iiave  tlu’  inverse  n-lationsliips 


G r 1’.'  ■ G, 


K = F.  • K, 


Carrying  out  the  operation  F • Cj  yiidds  the  inverse  of  Ec|s.(Z0) 


n 

iP  - V e,,  (fi[  , K.  • P'},  P-’  Z P' 

Ji  = !-'  1 ij  i-I-^ 


(23) 


(2-41 


Similar  results  are  obtained  from  li  • K,  with  L's  and  E's  la-placinL: 
ll's  and  P's,  rc'spectiv'cly . The  relationship  - ^ • G takes  the  lorm 


hi'- 

• 

•• 

. 7‘" 

z'^  • 

. z^'-” 

: 1 

q" 

-nl  .. 

"-j-nti 

>:-i  .. 

■ 

• 

I'iin  1 

7“  • 

■ 

l"l  ! 
1 

■ t 
ynn 

7'^'  ■ 

• w“‘7 

1 

■ 

L i*"  J 

where 


:ik 


ri 

= L' 

J - 


1 ‘ ij  Ki  O 


*^ik 


n 

= E 
J=i 


J-' 


>J  ' kj  t'.i 


R 


( ,h 
.1  k 


W 


ik 


11 

= E 
J=i 


ij 


C , ■ (H  , 
k 1 c .h, 
J 1 


'j-* 


R 


c .h 

J k 


W'iT) 


(23) 


(2b) 
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Situ'f  tho  rcterence  matrix  is  giMierally  sparser  than  the  path 
matrix  n , tiiere  arc  less  terms  involved  in  generating  the  elements  of 
V than  in  cimerating  the  i-lements  of  a . 

Transformed  Equations  of  Motion 

The  transformed  rnomentum  formulation  equation  of  motion  is 

G f X = K whtTe  X = A*^‘X-A  - G (27) 

and  the  I ransfornn.  d velocity  formulation  equation  of  motion  is 

! Y - K he  re  Y = A^  • ( Y+  a ■ A • ) (28 ) 

Performing  the  indiiaiti'd  operations  for  X and  Y yields 


X^  = 

y'  ^ 

: . v\  ■ 

• ■ V 

h • P" 

n 

d 

d j 

(29) 

n 

"^h 

- 

L (e;.  t-  R 

c.h . 

(30) 

j 

tfj 

1 

1 J 

a 

^-^c 

(31) 

t . 
1 

k f i i 

n 

c 

- L 

^c-  -h 

(32) 

if] 

k=i 

kj\ 

From  E',r|.  (H)^  for  v*,--  \vr  note  that  a^.-  is  the  p^t  of  the  acceleration 
v^,  which  is  not  linear  in  and  t^J.  r^vte  that  Ci  is  a ''Coriolis-type”, 

''i.i'nt  rifueal -ty  pe”  fore  i-  on  Body  i,  and  C-^  is  this  force  on  System  j. 

E^^  is  a co:  ibined  'Euler  coupling”  torque  j^us  monient  of  Coriolis- 
typi’  and  c ent  ril  ugal  - 1 ype  lorcc*.  Of_course,  X and  Y are  only 
fictitious  forces”;  the  true  force  is  K. 

Jerkovsky  (1977b)  showed  that  there  exists  a D such  that 

X -n  • Y ^ d‘  • “ (33) 

These  rel.tt 'onsbips  show  tlial  X and  Y are  closely  related.  However, 
n liivob.e-  . rv  ( omplicated  terms;  some  of  these  terms  drop  out  from 
f)  • , ,,nd  others  (Imp  out  from  t)*-  • Thus,  it  is  not  conceptually 

or  ( om[)utation-iIly  ell  ic  lent  to  cie  ter  mine'  D as  a prelude  to  dctc  r mini  ng 
X or  "Y;  i.  e.  , D is  not  a good  "intermediate”  or  "auxiliary"  variable. 

It  is  most  straightforward  to  cj^-termine  X from  Eq.  (29)^  and  Y from 
Eqs.  (29>2  to  (32);  i.  e.  , "vj^  , C*  , etc.  , are  good  intermediate  variables. 
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Tlu'  In  i-c|uaiions  in  Eq.  (Z7)  can  bn  written  out  as  follows 


(3-n 


Note  that  tlicsc  fciuations  arc  I'sscntiaUy  Eejs.  (1.6  1)  and  (1.  59)  of 
Mcirovitch  (1970),  written  for  System  i with  respect  to  the  :i  oviny 
point  h^.  Similiarly,  for  the  2n  eciuation^.  in  Etj.  (28)  we  havi- 


n rr 

V ( *f  J 
' h.h. 

1=1  1 j 

ii'  ■'  *S  , • 

CttIi. 

' f-i,  = 

} 1 

^ -t 

„h. 

■ • 

U-*)  a . P 

(35) 

j 1 B J 

The  tr.mstormed 

i-quitlions  of 

motion,  G t X 

= K ai'id  .a  • t Y - K, 

ri'  not  really  useful  in 

them  selves 

, because  they  < 

illow  6 degree's  of 

frei'dc^m  b<  twi-en  all  the  bodies  of  the  tree  c cjnfi  eu  r a 1 1 on.  If  tliere'  are 
really  t)  decri  es  of  freedom  between  all  the  bodies,  then  it  is  much 
sin, pier  to  use  the  priniiti\f  ‘-qxiations  of  motion  G + X = K or 

‘ Y K.  riie  re.'tl  utility  of  the  transformed  eejuations  is  that  in 
olitaininc  them  we  introdnii-d  a number  of  important  vi'Ctor  and  dyadic 
\.iriables  which  i .in  bo  used  ,i  s intirmc*diati'  or  auxiliary  variables  when 
the  final  ec|nalions  (with  U’ss  than  6 degrees  of  freedont  between  some  ill 
tlu>  Ijodies)  .ire  obtained. 


Separation  of  Free  and  Constrained  Motion 


To^r.ii  i I it.ite  tjie  imposUion  of  const  r.aints  in  relative  motion,  we 
exp.and  il'  , , llfij,  and  into  scalar  components  as  follows 


9! 

+ 

' 1 

y^2  ^2 

5 

/si 

T,i  i 

T ri 

u = 

u'i 

4 

2 

+ U, 

i 

1 

2_ 

3 

>'L 

.yP: 

^'h.,vP 

. y ^ '*  4 

n 

1 

1 

1 

i’  2 

“i 

4 

i 

1 

2 

3 

-^i 

6t 


(36) 


a, 

lor 

and 

yi- 

^ vV 


\ nd 


s 
1 s 


E^ 

to 

lor 


tlu 


• i ♦ 

are  c'xpanded  sim.ilarly  to  Hj.).  and  P , respectively. 
i arc'  3 gene;  rally  non-orthogonal  unit  vector  at  hj  in  Body  i, 
s - 1 to  3 are  3 orthogonal  unit  vectors  at  bj^  in  Body  i_. 
rocal  vc'ctor  to  , and  hence  7^,  ' • = 1 , t 


r ( ■ 1 1 p r 
= 0 if  r ^ s. 


g , ^s  ^s  **  and 

Now  introduce  the  following  column  matrices 
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rr 


\ 


— 

‘1 

i 

, A = 

*2 

y 

^3 

3 

— . — 

1 

■’ll 

II 

^2 

• ^A^  ^ 

^^4 

• «h.,r-''== 
1 

/k . 

(37) 


(38) 


Eqs.(36)  now  become 

, u' 

Now  making  use  of 


• t . 


. l 


T . , t r 

Ai’^h. 

= ‘‘'Psi 

(3b) 

• t 

= V a‘-  4'  =i3.r‘. 

(40) 

matrix,  we  can  invert  Eqs.  (39)  to  yield 

o'  - 

■ Pk  = ■ P‘ 

(41) 

= r^'"  • iV  , u i = fs 

' i'  1 

These  equations  can  be  inverted  again  , to  get  back  to  Eqs.  (39)  by  making 
use  of 

r^V'%T,  a^^a^=T,  (42) 


where  1 is  the  identity  dyadic. 

We  now  suppose  that  the  "gimbal  axes"  or  "Euler  angle  axes" 

9^  are  so  arranged  that  if  there  is  one  free  rotational. degree  of 
freedom  for  Body  i,  this  degree  of  freedom  is  about  7 J ; if  there  are 
two  free  rotati^.nal  degrees  of  freedom,  these  degrees  of  freedom  are 
about  9|  and  ■ Similarly,  we  suppose  that  the  "displacement  axes" 

0 g are  so  arranged  that  if  there  is  one  free  translational  degree  of 
freedom  for  Body  i,  this  degree  of  freedom  is  along  ; if  there  are 
two  free  translational  degrees  of  freedom,  these  degrees  of  ffeedom  are 

along  and  9^  . .This  orderly  separation  of  free  and  constrained  axes 

allows  us  to  write  , ri'”  , and  as  follows: 


r"  = 

1 

1 

T-.  i __ ! 
» ^ ■"  1 

i 

. A = 

i 

^f 

[ 

•-1 

o 

L 

1 

A^ 

c 

where  the  subscript  f denotes  the  free  axes,  and  the  subscript  c denotes 
the  constrained  axes.  We  can  similarly  separate  the  components  of 


A 
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, U . Hi,  , <uid  P 
i 


r - "1 

- . ” 

>4i’ 

“n , 

! 

_ 

, = 
A 1 

u‘ 

c 

1 

1 

Pj 

c 

Eqs.H-16)  can  now  be  expressed  as  (compare  with  Eqs.  (39)) 

= I 

f f c c 


(.•  C 


h!  = r;  'Si;  , 

n.  f i c c 
_ 1 

p' = ^|.’'p|.  t 

11  c c 


Inversi'ly  (compare  with  Eqs.  (41)) 

u|  “ r|-  ■ 


Q, 


U 


»;  - K ■ "i  ■ 
1 

,J  i 


h'  =r' 

c c 


u - f 

and  are  expanded  similarly. 


1^1  1 

P = A 
c c 


H 


h. 

1 


(44) 


(45) 


(46) 


V (.‘locity  T ransiormation 

/s  ^ /s 

We  are  now  in  a position  to  define^the  free  variables  '£  , Gf,  K£ 
and  thc'  constrained  variables  K^,  . We  define  A and  ■ as  follows 


^ /V 

' = A * 


(47) 


where  A is  a rectanijular  matrix  of  vectors  and  ' is  a column  matrix  of 
scalars.  Tiie  expanded  fortn  of  this  equation  is  as  follows  (compare 
v.ith  Eqs.  (45)) 
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•t 


(48) 


Note  that  as  indicatt'd  with  dashed  lines  in  Ec[.  (48),  this  ran  readily  be 
written  iti  partitioned  form  as  follows 


L A..  A 
t c 


/\ 

C 


/\  /\  /\ 

f ^ ■‘\- 


(49) 


The  inverse  of  this  relationship  is 


/s 

= n • 


B 

c 


1 

I f ; 


(50) 


when  B is  a rectangular  matrix  of  vectors.  The  expandi-d  form  of  this 
equation  is  ns  follows  (cornpan-  with  Eqs.  (46)) 
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1 Fn*  1 


[O',  idontly , A is  a matrix  of  unit  v<'Ctors.  B can  be  formed  by 
taking  the  transpose  of  this  matrix,  and  then  replacmg  c^ch  of  the  unit 
xect^rs  by  their  reciprocal  vectors:  symbolically^^  --  Note  that 

• A : the  bn  X 6n  identity  matrix,  whereas  aS  is  an  identity  with 

2n  1 's  on  the  diagonal. 

Induced  Transformations 


We  now  define  G and  K so  that  the  kinetic  energy  T and  its  time 
derivative  T maintain  the  general  form 


4 g'  • ^ = 


I pi  P 
" 2^ 


1 ^ ^ 
2-(Cf  ^ i G^^  ^ ) 


T = • 


This  requires  that  G be  defined  by 


/S  /S  f 

G = A^ 


r /N 

G = iG, 
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Carrying  out  the  indicated  operations  yields 
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c 

(54) 


Simila  r ccjual  ion  are  obt.iined  fc^  K,  with  L's  and  F's  replacing  H's 
and  P's.  Note  that  , Ci , and  K are  eolutiin  matrices  of  scalars. 


It  is  also  necessary  to  define  M-  by 
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(56) 


where  r,  s = f,  c,  and  where 


fj  . .TJ,  • r-'' 

r s r h.h.  s 

i J 

S^J  . s ^ ^ j' 

r s r c . n . s 

ij  1 

M^j  ^ . /i-i 

r s r s 

^ Z'  /V  ^ ^/\ 

TJbe  relationship  G = has  the  inverse  relationship  ' = •'G, 

where  ^ is  defined  by 


(57) 


/\  /\  — /Nf 

V =:  B . V . b'  = 


/s 

B 


L B,  B ] 
f c 


A 

«f  • 

nt 

^ 

«f  • 

V • 

-'^t' 

B 

t 

A 

'ff 

A 

fc 

A 

b 

c 

V . B^. 

A 

B 

c 

" 

'^t 

B 

c 

A 

cf 

A 

cc 

“ 

- 

- 

- 

Carrying  out  the  indicated  operations  yiidds 
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(591 


(60) 


Separated  Equation  of  Motion 

's 

N. 

'■V  The  momentum  formulation  equation  of  motion  is  now  given  by 

''  /S  /\  -A  A A f At  

G + X = K where  X - A • X - A • G (61) 

and  the  ve locity' *’orm ulation  equation  of  motion  is 
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ni’+ 


M-  -t  Y = K where  Y = A • (Y  + ►!■  • A (r. 

/\ 

It  is  actually  more  convenient  to  determine  Y from  Eq.  (f>2)2  with 
u ex['>anded.  Thus, 

Y = A • { A^  • L Y + hi  • (A  • ' + A • A ' ) j} 

Parviuoning  Eqs . (61)  and  (62)  into  free  and  constrained  parts,  yields 

/X  /\  /\  /Sf  — /\f  — 

G + X = K where  X = A • X - A • G (f 

r r r r r r 

z'-  ^ /V  A ^ '^t  ~ - ^ ^ 

u.  -thi.  ' +Y  = K where  Y A • (Y  + P • A - ) (6 

rl  f rc  c r r r r ' 

/N  ^ 

'•'.here  r f or  c,  and  where  X and  Y are  given  by 

r j.  a / 


. V • pi  - ■ pi 

r hj  r h 


• - 1'“'  • h" 


rl  • e‘ 

r h. 


X = 
r 


1 

A.  • C 


- a“  • 

r 


\ c 


whe  rt 


n 

= E 

i=I 

( E^  + R , ■ 

' c . c . h . 

c'  ) 

j 

1 1 J 

”►  • 

e' 

C 

= E^ 
c. 

+ 

c . 

'1  = i^k 

k € i 


k ei  i k 


= a 


■I^Note  that  A = 0 ; it  is  retained  in  X^  in  Eq.  (66)  merely  to  show  the 
general  pattern. 


(71) 


CJ  = L C‘  (71) 

ifj 

Note  that  M*  times  the  sum  ov^r  k for  the  first  term  in  Eq.  (70)  yields 
^ of  Eq.  (31)..  Also  note  that  5ri  is  the  p.a.rt  of  the  translational 
jet  deration  which  is  not  linear  in  and  • Similarly 

is  the  part  of  rotational  acceleration  which  is  not  linear  in  i2j-k 

and  U,|k  • is  clear  that  Xj.  is  conceptually  simpler  in  Y^. 

Eqs.  (64)  and  (65)  each  represent  two  equations:  one  for  the  free 

\ariables  (for  r = f),  and  one  for  the  constrained  variables.  The  constrained 
variables  can  be  dete rn'iined  algebraically  from  the  free  variables,  and 
therefore  in  a dynamics  simulation  we  only  need  to  "integrate  the  free 
variables  equation.  In  the  momentum  formulation,  the  required  differen- 
tial equation  is 


/V  /\  /N 

G^-  + Xj.  = 


touetlu  r with  the  algebraic  equation 


(72) 


(73) 


In  offi  ct,  il  is  necessary  to  invert  Uj-|-  ordc] 

the  "known"  quantities  Gf  (known  trom  "integration  ')  and  (prescribed). 
In  the  \elocity  {ormulation,  the  required  differential  equation  is 


n order  to  obtain  in  terms  of 
f/v 


/V  /N  /V  O ^ 

ff  f fc  c f 1 


(74) 


Note  that  in  this  formulation  w^e^nuist  also  effectively  invert  j...; 
in  addition,  we  must  generate  . 


The  constraint  force  K(.  i:^not  required  in  either  formulation. 
:1  it  be  desir 
lormulation  from 


Should  it  be  desired  to  uenerate  1^^,  it  can  be  done  in  the  momentum 


^ /S 

Ke  = G,  4 X^ 


,ind  in  the  velocity  formulation  from 

/\  ^ 


K = a’,  + 4.  4Y" 

( cf  f cc  c c 


(75) 


(76) 


.Note'  that  in  a vo-locity  formulatimr  simulation,  all  the  quantities  required 
to  gc-nerati'  the  constraint  force  are  already  "available  , whereas  in 
a momentum  formulation  it  is^iecessary  to  generate  G^..  We  will  give 
an  alternative  expression  for  Kf.  in  the  momentum  formulation  in  the 
next  section  (see  Eq.  (105)). 
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CiJniputational  Asp<  i t s 


III  lioth  till’  iiioiiu-iil  u!i  .iiid  voloiily  formulation  it  is  tioii-s-sa  r>  to 
fficctively  imurt  Ot  lourso,  ai  tual  in\iTsion  is  not  n-ally  lu-ii’s- 

sary  siiKu  it  is  only  ni  ii-ssary  to  "sol\i-" 


^ /N  ^ VS  /N 

— 

' (i  I ~ t ' ic  t 


(77) 


for  , or  to  "solvo” 


for  . Numi' r i (<i  1 ly  , thi-si*  two  "solving"  procnsses  are  identical 
and  some  sort  of  iteration  process  can  be  used  whic  h manipulates  the 
elements  of  However,  Russell  (1971)  has  sho\.^n  that  it  is  more 

effii  ient^to  compuli^otily  the  "diagonal  blocks"  of  a and  to  put  the 
i-est  of  . tine  s j-  on  the  right  hand  side,  without  explicitly  coniputing 
the  rest  of^the  elements  of  This  form  of  "block  iteration"  is 

discusser!  by  V.ii’ga  (TtbiZ). 

Another  method  of  ai  oiding  the  inverse  ot  a j-|-  is  to  compute  the 
inverse  of  ^ and  then  use  the  relationship 


Cv  Cs  /N  ^ 

‘Nf  f-^'f-  Yf-"fc  'c 


(78) 


/N  - 1 /S 

"ff  " "ff 


/\  /N.  J /X 

w V V 

f c c c c f 


(79) 


Ignorint;  the  constraints,  there  arc-  6n  degrees  of  freedom  in  a system  of  n 
rigid  bodies.  If  nf  of  the  degrees  of  freedom  are  free  and  n^-  are 
constrained  ( iij-  t n^,  --  6n  ) then  the  use  of  Eq.  (79)  may  be  i^sirable  if 
n^.  nj-  , or  it  might  lx-  desirable  to  use  Eq.  (79)  because  may  be 
sparser  (or  otherwhsc'  simph-r  to  compute  because  of  less  multiplications) 
than  . 


Another  aiternafi\i>  to  analytically  gent;rating  all  the  elements  of 
a i s to  use  the  ri’lat  ionship  ^ in  the  form 

G • _ a • ( ^4  )j  where  = A • A (80) 

The  j**^column  of  a is  then  obtained  by  letting  the  clement  of 
in  Eq.  (80)  be  unity  ,^while  the  rest  of  the  elements  of  are  zero, 
and  then  computing  G in  thn-e  stages  as  indicated  via  parenthesis  and 
bracket  in  Eq.  (80).  The  6 obtained  in  this  manner  is  then  the 
column  of  "a  . This  aiiproach  has  been  used  by  Russell  (1969).  The 
elements  of  - can  be  obtained  similarly  from  the  three-stage  computa- 
tion 

R ■ [ '-  • G ) j where  R = B'  B (81) 

X ^ -1 

and,  of  course,  R = x . Some  computation  and  computer  storage  space 
c an  be  saved  because  g"  and  are  symmetric. 
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Strictly  speakirijj,  the  equations  of  Russell  differ  somewhat  from 
the  mon.entum  formulation  equations  presented  herein,  because  Russell 
uses  velocities  and  momenta  relative  to  the  composite  center  of  mass. 
Such  equations  can  be  obtained  by  transforming  from  the  abo_^e  trans- 
formed velocity  ' to  “ which  is  the  same  as  “ except  that  is 
replaced  by  Vf-j  , the  velocity  of  the  composite  c^ter  of  mass.  " ^ 
then  replaced  by  a similar  ; similarly  G'and^^'are  replaced  by  “U  ' 
and  K . In  fact,  Russell  (1969)  refers  to  G'  and  G ' as  "primed  momen - 
turn".  In  such  a formulation  one  also  obtains  a M.' , ' , etc.  , but  for 

computational  efficiency,  Russell  always  avoids  analytic  generation  of 


Whether  a momentum  or  velocity  formulation  is  used,  the 
"integration  accuracy"  must  be  checked.  One  time-honored  check  is 
checking  the  constancy  of  the  inertial  i omponents  of  total  angular 
momentum  during  periods  of  zero  external  torcjues  A m^e  comprehen- 
sive check  would  be  to  compute  and  integrate  T = ^ 

compare  this  integ rated  value  of  kinetic  energy  with  the  kinetic  energy 
obtained  J^^im  T = y (uj-  f G^^c)^  Note  that makes  a contribution 

to  T if  ' is  not  zero;  similarly,  G contributes  to  T. 

c 

Coupling  of  Free  and  Constrained  Motion 


As  an  alternative  to  transforming  to  free  and  constrained  variables 
we  can  leave  the  equations  of  motion  and  equations  of  constraints  in  a 
coupled  form,  and  then  solve  the  equations  of  motion  and  equations  of 
constraints  simultaneously.  Recall  we  started  out  with  G + X = K or 
a • ■ ^ Y = K and  then  inade  the  transformation 


/\/s 

J 0 


whe 


J = A ■ A 


(82) 


This  can  be  written  as 


f 


' 

f c c 


(83) 


wht 


/X 

4 = A 
s 


■ A for  s = f and 
s 

' = R • - 


c.  Inversely,  we  have 


/N 

Rr- 

. . _ 

f 

f 

/N 

Rc| 

(84) 


'ts  ^ 

where  R = B 


B and 


/\ 

R 


XS 

B 


B for  s = f and  c.  Since  we  consider 
to  be  prescribed,  we  effecuvely  have  the  equation  of  constraint 


/N 

R 


/N 

c 


(85) 


There  are  n^  scalar  elements  in  and  thus  th«-re  are  n^,  scalar 
equations  of  constraints. 
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From  the  equation  a - A (t  we  get  the  equation 

/N  /V 

K = J • K = 


1 

II 

1 

I 

1 

1 

1 

n 

l 

Inversely,  we  have 
K = R K = 


Thus,  we  can  write 


^t  /Nt 


/N,  ^ /\i  /N 


(86) 


(87) 


K = 


where 


K"  = K, 


= R K 
c c 


(88) 

(89) 


is  the  ''applied''  part  of  K,  and  is  the  "constraint"  part  of  K. 


If  the  multi-body  configuration  is  not  ^tree,  it  is  not  a simple 
matter  to  find  the  transformation  operator  A and  the  appropriate 
velocity  ^ with  free  elements  and  constrained  elements  ' . 

However,  it  is  usually  simple  to  obtain  an  expression  of  the  form 


c 


c 


(90) 


for  the  constraints  of  the  multi-body  configuration.  Here  R^,  is  some 
(primitive)  velocity  constraint  transformation  operator,  and  is  some 

prescribed  constraint  velocity.  In  a tree  configuration,  we  simply  take 
to  be  and  we  take  ' to  be  . From  the  general  properties 

of  Lagrange  multipliers,  it  now  follows  that  we  can  write 


K = where  K (91) 

c c 

where  is  the  Lagrange  multiplier  column  matrix_^(usually  denoted 
by  ^ ).  In  the  case  of  a tree  configuration  we  take  to  be  K^,.  Thus, 
the  momentum  formulation  equations  take  the  form 

G + X = + R^';  (92) 


and  the  velocity  formulation  equations  take  the  form 

u • • + Y = + R^  K 

c c 


(93) 


In  either  formulation  we  also  need  the  constraint  relationship  in^ 

Eq.  (90).  We  will  now  examine  what  is  required  to  determine  K^-  in 
either  formulation. 
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Velocity  Formulation 


In  the  velocity  formuiation,  the  simultaneous  equations  which 
must  be  solved  are 


c 

k''  -y 

R o 

c 

. V • 

(94) 


where  the  second  row  of  this  matrix  equation  was  obtained  by  taking 
the  time  derivative  of  Eq.  (90)  . Inverting  this  relationship  yields 


! 

• 

j a b 

- K 

! c d 

wh  e r e 

1 

d - - ( R^,  ■ 

' K“  - Y 


. I 


0,  • 


t ,-l 


a = '“  + ^ • 

b = ~ • 

c = ' d R 


R ) 


d R 
c c 


R'  d 

t 


(95) 


(96) 


1 i e n c e 


K = - c 
c 


(K‘ 


Y ) - d ( - R^_ 


The  differential  equation  of  motion  is  now 


(97) 


==  a 


( K“  - Y ) + b ( - R^ 


■ ) 


Note  that  for  a tree  configuration  we  may  take-  R 
1 is  equal  to-'-'J,  . The  evaluation  of  -'3.  tan  of  *i- 


then  d is  equal 
by  an  evaluation  of  according  to  the  relationship 

^ ^-1 


1 /X 

- U 

c c c c 


cf  '^T  fc 


(98) 

to  be  R , and 
ourse  be  replaced 

(99) 


Momentum  Formulation 


In  the  momentum  formulation,  the  simultaneous  equations  whuh 
must  be  solved  are 


t 

» f^c 

~ 

r 

- 

G 

X 

1 

6^-  ^ O 

- K 

c 

c - fie' 

- R,--  G 

- 

- 

(100) 
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where  the  second  row  of  this  matrix  equation  was  obtained  by  setting 
= V • G in  Eq.  (90)  and  then  taking  the  time  derivative.  Inverting 
this  relationship  yields 


G 

i ^ 

b 

K'"'  - X 

• 

• ^ 

-K 

i c 

(i 

- R , ■ 

- R,  • 

V • G 

c 

c c 

C 

whe  re 


d = - ( PI 


B*  - d 
e 


b = 


I r d R 
c c 

(V  d = a 


V = 4 


c - - d 


V = c 


(101) 


(102) 


Hence 


K = - 
c 


(K^  - X)  - d (~  - R^ 


- R. 


(103) 


The  equation  for  G in  Eq.  (101)  is  actually  less  convcmiimt,  though 
equivalent,  than  the  original  i*quation  G = - X + R^  K^-.  The 

equivalence  of  Kj,  in  Eqs.  (97)  and  (103)  follows  from 

Y = X+  '.I  ■ =X-u-v.G 


(104) 


Thus,  whether  a velocity  or  a momentum  formulation  is  used,  essentially 
the  sanu-  ecjuations  ar^  involved  in  determining  the  constraint  force  (or 
Lagrange  multiplier)  K^. 

/s  . , 

Recall  that  when  we  discussed  solving  for  K_  in  the  previous 
section,  we  found  that  Eq.  (75).  for  the  momentum  formulation  was 
not  really  convenient  because  6^,  is  not  available.  We  now  sec  that  we 


(105) 


can  use 

11 

^ - 1 ^ 

- V 1 R . 

cc  t 

' - 

■ <^f  Kf 

^-1  A 

- Y)  + 'e,  ( 'e 

A 

- Rc- 

■ ) 

= 

/N_  1 yv 

"cc 

- • 

/St  /s 

( Rf  Kf 

- X)  + (^ 

' CC  c 

A 

- Re' 

-Rc  • 

V • G) 

Note  that  either  forms  of  Eq.  (105)  could  be  used  wdth  either  the 
velocity  or  the  momentum  formulation. 

Comparison  With  The  Literature 

Table  I lists  some  of  the  better  knowm  references  on  multi-body 
spacecraft  dynamics  and  briefly  comments  on  them  in  the  light  of  the 
above  discussion. 
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'Dll'  i (|uatjons  of  Hookor  and  Margulies  (1965)  and  of  Roberson 
and  \\iiti  nburL;  (Idno)  (.an  be  obtained  in  two  stages  as  follows:  First, 

i e.\|jress<.‘d  in  ti-rms  of  _ which  consists  of  the  inertial  angular 
velocities  of  each  body  plus  the  inertial  translational  velocity  of  the 
coir.po.site  centi>r  of  mass;  the  i‘qi»ation  of  motion  then  is 
••  • T V - K where  - A|  . ' , a --  A*  • • Aj  , and  K ~ A]  • K;  the 
elenu'iits  (,)f  ^ can  lu'  expressed  in  ti'rms  of  " bar y c ent c r s"  and 
augnu-nti-d  bodies".  Second,  the  Iransfor tru‘d  equation  of  motion  is 
' oupUai  With  tlu'  m'lativa-  rotation  c onstr.iint  ecjuation,  B2  • =0,  and 

the  constraint  torcjue.s  are  obtained  \ia  Lagrange  multiplic.' rs . 

I o relate  the  refereni  (.-  matrix  r to  the  "incidence"  matrix  S of 
Kohc'rson  and  Wittenburg,  let  ■ | b(.>  the  first  row  of  r , and  let  C'2  be  the 
remaining  n-l  rows.  Ihen  - S* , and  eacli  column  (.if  S (or  row  of 
2)  (.ontains  all  zero  elimients  except  for  oni'  t 1 and  one  -1.  Similarly, 
tne  jiath  matrix  can  be  la-lated  lu  the  "ordering"  matrix  T of 
Roberson  and  Wittenburg  (somelinu's  denoted  by  S-=)-  Let  be  the 

lirst  ( iilumn  of  " , and  K‘1  ^ bi’  tlie  remaining  n-l  columns.  Then 

2 T*  , and  the  lirst  column  ol  T (or  row  of  ’'2)  contains  all  zero 
ehnnents.  .Sini'e  and  - are  inv(.'rsi’S  of  each  other  ( c-  1^,  the 
n ■ n ideiiiity  matrix),  it  follow.s  that  ^ ”2  ~ *n-L 

inverse  of  consi’ciueni  ly , .S’-  d’  - ln-l>  on  TS  - In-L  ^ 

'I'he  approach  of  Velman  (1967)  is  similar  to  that  of  Mooker- 
-Margulies  and  Robe  rson- Wittenburg,  e.xcept  that  inertial  angular 
veloc'itii's  are  replaced  by  relative  anguLtr  velocities  (and  relative 
liiietir  veloiities  in  case  of  point  masses).  The  equation  of  motion, 
befyre  ijnposition  of  the  relative'  riitational  constraint,  has  the  form 
- K , when'  , 'S'  , and  K are  column  m.itriecs  of  scalars, 
j^nd  ci  is  a (positive  delinite  synmu'tric)  matrix  of  scalars.  Next, 
and  its  inva-rse  B are  introduced  as  follows:  A and  _ - B 

P;i  rtitioning  A and  B (similarly  to  the  [lartitioning  of  A and  Pi  in 
Eqs.^(49)  and  J50))^qelds  Aj-  'f  t A,.  ' and  ‘|-  Bj-(r.  A ^ 

K = Aj^K^and  B’,K  are^simi  !a  rly  pa^rtitioned  as  K|  K > 

Kj.  “ a’.  anjl  K = B|;  K t B|^,  K - K‘^  ; K^-',  where_K‘^  is  the  "ajiplied" 
part  of  Kand  is  the  "lonst  raint"  part.  From  AB  - I there  follows  thi 

relationship  f ; 1,  whc'rg'  1 is  a bn  'jfin  scalar  identity  matrix, 

and  whi're  'q  - > ' e R A = 1 tht'ri'  follows  thi' 

relationships  ^fAj-  - Ij-,  -■  1^.,  BfAg  and  l\  A]  O^.j-  : 

here,  1^-  and  1^  are  scalar  identity  matrici's,  and  0_f^  and  0^,|-  an- 
matrices  of  zeros  (which  are  transposins  of  each  other).  and  are 
6n  ■ 6n  matrices  of  scalars  which  are  idempotent:  'i^-  ” ;ind 

'■(.  '!(-  - *’q-  ; hence,  these  matrices  are  "projector. s ' (or  "projection 
operators").  In  the  case  of  Velman  (1967),  A nnd  B are  permu'.-ition 
matrices  and  therefore  B = A*-;  consequently,  Bj-  = Aj-  an(J  B^^.  - , 

and  it  follows  that  ajid  ^ are  symnu'tric.  ^ In  terms  of  these  pre- 
je(  tors,  the  ecjuation  Aj.  = 0 bmairnes  fj.  K*'  = 0 upon_^left  multi- 
plying by  b|.  , and  the  g'qua tion  q,  = B^,  ~ becomes  ^ = '^^c~c 

left  multiplication  by  A^. . Now  that  Velman' s symmetri(.  projei  tors 
have  been  introduced  in  terms  of  the  transformation  operator  formali.sm 
the  rest  of  Velman' s procedure  can  be  set'll  from  thi'  dis(  ussion  by 
Likins  (1970). 
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ative  translation; 


1 hi’  I um  pul.  1 1 1 ona  1 .lsp^lt^l  ol  llu  1 1 ooku  r - Mu  r g ul  i c‘ s and  Hoberson- 
\V  itti'nburg  formalisms  .iro  tri'uti’d  in  ]•' Ui  si'lu' r (1971)  and  tarrell  et  al. 
(19t)8),  ri’spi’i  tivi’ly  ■ l-'l  I’i  si  ho  r al  so  di  s c us  sos  Volman's  procedure 
for  I’liminatinu  tlio  i onst  r.ii  nt  s . 

The  equations  oi  Palmer  (19o7)  were  the  first  general  set  of 
equ.itions  wlu're  the  l onstraint  torques  .ire  di'coupled  by  transformation, 
ratlu'r  than  solved  via  Lagrange  multipliers.  However,  Palmer's 
equations  are  restricli’d  to  configurations  which  are  "clusters"  in  the 
sense  that  all  of  the  bodies,  other  than  Body  1,  are  attached  to  Body  1. 

It  IS  interesting  to  note  that  Palmer  (1967)  discusses  velocity,  accelera- 
tion, and  torcjue  transformations,  but  he  does  not  make  extensive  use 
of  these  t r.i nsfo r ma t i ons . 

Russell  (1909)  was  the  first  to  develop  a set  of  transformed 
equations  for  an  arbitrary  tree  configuration.  Russell  chose  a momen- 
tum formulation  where  the  dynamics  state  variables  are  the  free  compo- 
m-nts  of  the  transformed  momentum;  however,  the  constrainc-d  compo- 
nents of  momentum  and  the  primitive  and  transformed  velocities  are 
retaini'd  ;i  s i nt  e rmediat  t'  or  ".auxiliary"  variables  so  that  the  final 
equations  have  a particularly  simple  form. 

Farrenkopf  (1969)  introduced  an  "inductive"  method  of  "digitally 
synthesizing"  the  dynamics  equations  via  a "c  ombining  algorithm".  The 
original  formulation  was  restricti'd  to  tree  configurations,  but  it  was 
extended  to  ti'rminal  flexible  bodies  by  Ness  (1971).  Their  equations  are 
essentiall>  the  "transformed"  equations  of  this  paper  if  left-multiplied 
by  a non-singular  matrix;  this  left-multiplication  is  necessary  to  make 
the  Farnnko[3f  "mass  m.itri.x"  synimetric,  as  it  is  in  this  paper.  Thus, 
if  P'a  r renkopf ' s equation  of  motion  is  in  ' + y = k , then  left-multiplying 
6y.  yields  u ' f ■^  = K where  j m , Y=  »^y  , and  K = k. 

If  it  is  the  inverse  of  fv  , then  k = a*  K ; thus,  c.an  be  identified  as 
the  matrix  which  tr.ansforms  the  transformed  force  K of  this  paper  to 
the  force  k of  Farn-nkopf;  the  matrix  which  symmetrizes  m by 
left-multiplication  is  then  the  inverse  of  this 

The  Hooke  r-Margulies  formalism  was  converted  to  an  approach 
which  uses  relative  gimbal.angle  rates  in  Hooker  (1970).  The  resulting 
equation  of  motion  is  a'  * ' + Y'  = K',  where  K'  is  the  "primed  force" 
introduced  earlier  in  connection  with  Russell's  "primed  momentum" 
approach.  It  is  interesting  to  note  that  in  the  Hooker  (1970)  formalism, 
Russell's  primed  momentum  is  simply  G'  - p ' ' ' , the  kinetic  energy 
is  T - G'*^  It'  , and  the  time  derivative  of  kinetic  energy 

is  t = k't  *'  . ^ 

Likins  (1973)  extended  the  Hooker  (1970)  equations  by  allowing 
the  terminal  bodies  to  be  flexible,  and  by  allowing  each  rigid  body 
to  contain  axisymmetric  rotors. 

The  Roberson- Wittenburg  formalism  was  converted  to  an  approach 
which  uses  relative  gimbal  angle  rates  in  Roberson  (1972)  and  Wittenburg 
(1973).  Both  of  these  extensions  allow  relative  translation  between 
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bodirs,  tlu‘  extension  by  i^obersoii  also  allows  the  bodies  to  be  deform- 
able. but  the  equations  ol  motion  for  the  deformation  coordinates  are 
not  pr('sented.  Boland,  Satiiin,  and  Willems  (1974,  1975)  have  also 
obtained  Roberson- Wittenburq  ty[)e  equations  in  relative  gimbal  angle 
rates  and  relative  translational  rates,  and  have  described  the  use  of  this 
formalism  in  configurations  with  c losed  loops  (by  cutting  as  many  loops 
as  required  to  form  a trc'e,  and  then  introducing  constraints  via 
Lagrange  multipliers). 

The  Hooker  (1970)  eejuations  have  been  converted  to  a set  of 
equations  not  using  "barycenters"  and  ''augmented  bodies"  by  Frisch 
(1971,  1975),  Ho  (1974,  1977),  and  Hooker  (1975).  In  this  latest  version 
the  dynamics  state  variables  are  the  relative  gimbal  angle  rates  between 
bodies  plus  the  inertial  linear  velocity  of  a material  point  of  Body  1.  The 
resulting  equation  is  the  equation  ii  a + Y = K of  this  paper  (except  that 
these  papers  by  Frisch,  Ho,  and  Hooker  do  not  allow  relative  transla- 
tion between  bodies  - -howeve r,  the  terminal  bodies  are  allowed  to  be 
fle.xible;  Frisch  (1975)  treats  all  bodies  as  flexible). 

The  approach  of  Ho  (1974,  1977)  and  Hooker  (1975)  has  been 
e.xtended  to  a chain  of  flexible  bodies  in  Ho,  Hooker,  Margulies,  and 
Winarske  (1974).  The  most  interesting  feature  in  this  extension  is  the 
use  of  quasi-static  modes  plus  vibration  modes  to  describe  the  deforma- 
tion of  the  flexible  bodies;  the  use  of  these  modes  allows  decoupling  of 
the  constraint  forces  and  torques. 

The  formulation  of  Bodley,  Devers,  and  Park  (1975)  is  the  most 
general  of  those  in  Table  I.  It  allows  all  bodies  to  be  flexible,  it  allows 
up  to  six  degrees  of  freedom  between  bodies  and  any  of  these  degrees 
oi  freedom  may  be  prescribed  function  of  time,  and  it  allows  closed 
loops.  The  dynamics  equations  are  retained  in  "primitive"  or  "free 
body"  form,  and  the  constraint  forces  and  torques  are  obtained  via 
Lagrange  multipliers.  It  is  interesting  to  note  that  Bodley  et  al.  (1975) 
make  fairly  explicit  use  of  velocity  transformations. 

Of  all  the  authors  in  Table  I,  only  Russell  uses  a momentum 
formulation.  The  transformation  operator  formalism  was  initially 
developed  in  terms  of  a momentum  formulation  (Jerkovsky,  1976),  and 
the  extension  to  a velocity  formulation  was  made  in  order  to  provide 
an  overview  of  the  alternatives  in  Table  I.  As  a matter  of  record,  it 
can  denoted  that  the  use  of  a momentum  formulation  is  also  advocated 
by  Bodley  and  Park  (1972),  and  by  Williams  (1976). 

Conclusion 


An  overview  of  the  structure  of  several  multi-body  dynamics 
formulations  has  been  presented  in  the  language  of  the  transformation 
operator  formalism.  The  following  alternatives  have  been  discussed. 

1.  Momentum  or  velocity  formulation 
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2.  Si'p.i  iM li iii:  thi-  (.'qu-ilionti  ol'  motion  from  tho  equations  of  con- 
straints, or  i.  ouplin<4  tht-se  o(|uations 

5.  ComputniLt  tlu’  entire  mass  matrix,  or  computing  only  its 
l)lo(.k  cli,iL;C)nal  elements 

■1.  Compuline  tile  mass  matrix  analytically  or  numi'rically 

S.  Invertine  aj-j-  or  (both  ari‘  positive  definite  symmetric) 

It  has  bei  11  noted  that  the  same  ty|>e  of  linear  simultaneous 
ec|uations  must  be  soh'cd  in  the  momentum  formulation  and  in  the 
\eloiity  formulation.  The  same  mass  malri.x  and  the  Sc'ime  force 
apjjear  in  i-ither  formulation.  In  fact,  the  same  mass  matrix  and  force 
ari'  obtained  in  the"niatrix  method"  of  structural  analysis  if  the  trans- 
formation matri.x  is  i ontinuously  upd<iti*d  to  reflect  the  instantaneous 
\alues  of  the  coordin.ites . However,  tlii’  mat  rix^method  of  structural 
analysis  does  not  gene  rati'  thi-  extra  ti'rm  X or  Y bectiuse  the  time 
derivalivi'  of  the  transformation  matri.x  is  neglected. 
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